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ABSTRACT 

We challenge the conventional view that the majority of Mgll/Lyman limit absorbers 
are extended halos of galaxies comprised of "clouds" with near-unity covering factor. 
The gaseous disks of spiral galaxies are known to extend to large radii and likely con- 
tribute a significant cross-section for absorption. We perform a Monte-Carlo survey 
of QSO fields in which the model galaxies have Mgll absorbing "clouds" in a spheri- 
cal halo or in a randomly oriented disk. For both geometries, models that recover the 
observed properties of Mgll absorbers have only a 70-80% covering factor. Therefore, 
regardless of absorber geometry, a survey of randomly selected QSO fields should yield a 
non-negligible number of non-absorbing galaxies at small impact parameters. Very few 
have been observed (Steidel 1995). However, selection effects are important, and once 
the observational procedures are applied to our model fields we find that this result is 
expected. Since both spherical halo and disk models can be made consistent with survey 
results, we present tests for discerning the geometric distribution of absorbing gas from 
Hubble Space Telescope images and high resolution spectra. 

Subject headings: quasars: absorption lines — galaxies: structure — galaxies: evolution 
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1. Introduction 

Quasar absorption lines (QALs) provide unique, 
sensitive, and direet probes of the gaseous content in 
galaxies over look-back times comparable to a Hub- 
ble time. Studies of QALs promise to contribute in- 
sights into the processes that regulate the evolution of 
galaxies and the mechanisms that determine the phys- 
ical conditions in interstellar and halo clouds. The 
gas in galaxies is a key ingredient to their on-going 
evolution, yet virtually nothing is known about the 
gaseous conditions in early-epoch galaxies, and how 
these conditions evolve to the present time. The field 
of QSO absorption lines has matured quite rapidly 
in the last half-decade. As such, a reasonably sized 
database, from which we can draw inferences about 
the general properties of galaxies known to give rise 
to absorption, has accumulated. 

Low resolution spectroscopic surveys of the reso- 
nant CIV AA1548, 1550 and Mgll AA2976, 2803 dou- 
blets (Sargent, Boksenberg, & Steidel 1988, Sargent, 
Steidel, & Boksenberg 1988, Steidel & Sargent 1992) 
have established the statistical nature and evolution 
of various "populations" of absorbers. The evolu- 
tion of the population selected by the presence of 
Mgll is virtually indistinguishable from the popula- 
tion selected by the presence of a Lyman limit break. 
This comes as little surprise, since Mgll is known 
to trace HI gas that is optically thick at the Ly- 
man limit (Bergeron & Stasinska 1986, Steidel & Sar- 
gent 1992). Follow lip imaging surveys of Mgll ab- 
sorbers have been highly successful at detecting a 
galaxy at the redshift seen in absorption (Bergeron 
& Boisse 1991, Steidel 1995). The associated galax- 
ies are seldom fainter than 0.1 L|f and span the full 
range of galaxy types (Steidel, Dickinson, & Persson 
1994, Steidel 1995, Drinkwater, Webster, & Thomas 
1995). However, galaxies in the same cluster as the 
QSO provide an exception to the rule that every 
bright galaxy produces Mgll absorption (Bechtold & 
EUingson 1992). This could occur more generally in 
galaxy clusters where ram pressure or high energy ion- 
ization fields act to destroy extended low ionization 
gaseous components. 

It has become apparent that these various popu- 
lations are likely the absorption signatures of various 
regions and conditions within galaxies. Indeed, the 
survey data have been used to infer the statistical 
properties of a "typical" Mgll absorbing galaxy at 
intermediate redshifts (Steidel 1993a). Such a galaxy 



is characterized by an inner region of radius ^15 kpc 
that gives rise to damped Lya lines, a region ex- 
tending to ~ 40 kpc that produces Mgll absorption 
and a Lyman limit break, and an outer region ex- 
tending to ~ 70 kpc that is less shielded from the 
extragalactic background radiation and produces ab- 
sorption lines in higher ions such as CIV. Beyond 
these inferences drawn from statistical cross-section 
arguments, there has been little direct information on 
which parts of galaxies give rise to Mgll absorption. 
Studies through the Galactic Halo (Savage 1993, Sav- 
age, Lu, & Sembach 1995, Sembach, Savage, & Lu 
1995) have provided results which suggest that com- 
plex dynamical disk/halo interactions may give rise 
to Mgll absorbing gas. These processes may have 
analogs in external galaxies, but definitive interpreta- 
tions of absorption from the Galaxy are lacking due 
to our difficult vantage point. Bowen, Blades, & Pet- 
tini (1995a, b) examined UV spectra and images (of- 
ten in several bands, including X-ray and HI) of local 
Mgll absorbing galaxies. They often found that it is 
a challenge to unambiguously identify the substruc- 
tures giving rise to Mgll absorption in these nearby 
systems, which often showed some sign of past inter- 
action or were accompanied by satellite galaxies. In 
other words, the environments around the galaxies 
were not always isolated enough to clearly determine 
with which object the gas was actually associated. 
In the case of a line of sight through M81 and the 
Galaxy, each of which have satellite galaxies affect- 
ing their environments, the Mgll absorption profile 
is blended with a 400 km spread. A further im- 
portant lesson from their study is the possibility of 
altogether mis-identifying the absorbing galaxy. In 
the case of Q1543 -|- 489, two galaxies with veloc- 
ity separation 120 km s~^ are absorber candidates. 
The first has impact parameter 45^."'^ kpc, and the 
second 83/i~^ kpc. This latter galaxy is outside the 
range probed by existing intermediate redshift imag- 
ing surveys, and as such, would be completely missed. 
Indeed, such surveys would likely claim they had suc- 
cessfully identified the former galaxy as the absorber. 
Such mis -identifications or ambiguities, even if they 
are relatively infrequent, may have an impact on in- 
ferences drawn from small number statistics, such as 
the fraction of non-absorbing galax;ies at small impact 
parameters. 

Ultimately, the kinematic-rich high resolution ab- 
sorption spectra obtainable with 10 meter class tele- 
scopes and the high spatial resolution Hubble Space 
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Telescope (HST) images should aid the quest to un- 
derstand (in individual cases) which types and parts 
of galaxies give rise to Mgll absorption. Indeed, high 
resolution spectra have provided some leverage for ex- 
amining statistical correlations between host galaxy 
properties and absorbing gas kinematics (also chem- 
ical and ionization conditions). Lanzetta and Bowen 
(1992) studied five absorbing galaxies at resolutions 
of 7 and 35 km s~^ and found some profiles to be 
consistent with the signatures expected from a sim- 
ple rotating structure and some that match predic- 
tions from infall/outflow within a spherically sym- 
metric structure. Additionally, they found prelim- 
inary evidence that the profile complexity and the 
strength of absorption are correlated with impact pa- 
rameter (eg. only one or two subcomponents at larger 
impacts). However, as concluded by Bowen, Blades, 
& Pettini (1995b) from their study of nearby galaxies 
(where high spatial resolution data are available and 
evidence for interactions has been seen), generalized 
models of halos may be misleading. Detailed models 
of many individual systems may provide a database 
from which a more general understanding of absorb- 
ing gas can be obtained. Churchill, Vogt, & Stei- 
del (1995) presented a preliminary study of a striking 
HIRES (Vogt 1994) profile from Q1331 + 170, and 
suggest that quasi-symmetric outflowing structures 
(eg. supcrbubbles) may be present in what appears 
to be two close galaxies. Unfortunately, such a model 
can never be demonstrated to be a unique interpre- 
tation of the data. Clearly, even with HST images 
to corroborate inferences based upon high resolution 
spectral data, a large sample will need to be com- 
piled before a comprehensive understanding of the 
spatial distribution, motions, and origins of galactic 
substructures giving rise to Mgll absorption is devel- 
oped. 

Nonetheless, now that a database of Mgll absorb- 
ing galaxy properties (eg. luminosities, optical/IR col- 
ors, redshifts, and impact parameters) is becoming 
available (Steidel, Dickinson, & Persson 1996, here- 
after SDP), we can begin to learn something about 
the absorbing gas and its relation to the host galaxy. 
Though the high resolution spectra provide detailed 
information about the absorbing gas, existing low res- 
olution absorption line data are proving to be suffi- 
cient for developing a basic appreciation of the condi- 
tions under which the presence of Mgll is expected. 



1.1. This Paper 

In this paper, we are motivated by the basic ques- 
tion: "what is the general geometric cross-section of 
Mgll absorbing gas in and around galaxies?" What 
we really aim to learn is if a single geometric cross - 
section can be invoked to help predict the presence 
or non-presence of absorption once a galaxy's general 
properties and line of sight impact to the QSO are 
known. The conventional wisdom is that Mgll ab- 
sorption is produced in effectively spherical extended 
halos [dating from Bahcall & Spitzer (1969)]. Recent 
support for this view has been presented by Steidel 
(1995) (hereafter S95), who has demonstrated that 
the physical extent of Mgll absorbing gas around 
galaxies scales with rest K luminosity following a 
Holmberg-like relationship (Holmberg 1975) 

i?(LK) = 38/1-1 (i^AJcf j^p^ 

(cf. Fig. 2 in S95). The S95 paper presents prelimi- 
nary results from the extensive study of SDP which 
is still in progress. One striking result to date is 
that of 58 absorbing galaxies and 14 non absorbing 
galaxies measured, there are only two non absorbing 
galaxies "below the R{Lk) line" given by this rela- 
tionship (i.e. only two galaxies have line of sight im- 
pact parameters that fall within the now predicted 
surrounding absorbing region and do not exhibit ab- 
sorption). Also, in only three cases is the identified 
absorbing galaxy found "above the R{Lk) line". The 
most straight forward interpretation is that Mgll ab- 
sorbing gas is distributed with a roughly spherical 
cross-section, and has a nearly unity covering factor 
down to the Steidel & Sargent (1992) survey equiv- 
alent width limit lVo(A2796) ~ 0.3 A with a K lu- 
minosity (mass) dependent "cut-off boundary" . This 
conclusion is supported in part by the notion that a 
disk-like geometry, of random orientation, presents a 
significantly smaller statistical cross-section such that 
a non-negligible number of non-absorbers below the 
R{Lk) line would have been detected. 

However, these simple geometric models deserve 
further consideration, since there is clear evidence for 
"dumpiness" in the absorbing gas. Though Lanzetta 
& Bowen (1990) found that IVo(A2796) and the im- 
pact parameter D of the absorbing galaxies follow the 
"smooth" relation Wa oc D"" with a = 0.92 ± 0.16, 
S95 has shown that in fact there is a large scat- 
ter, which is not consistent with the tight correla- 
tion expected to arise from a smooth distribution of 
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halo gas (cf. Fig. 3 of S95). Since intermediate and 
high resolution absorption lines typically exhibit mul- 
tiple absorption subcomponents (Petitjean & Berg- 
eron 1990, Lanzotta & Bowcn 1992), we infer that a 
near-unity covering factor must result from absorbing 
gas distributed in discrete clumps which each have a 
minimum equivalent width Wmin set by the physics 
of photoionized diffuse gas. We also note that the ev- 
idence for unity covering factor is statistical - in any 
single case it is not possible to be certain that the 
galaxy identified at the same redshift as that seen in 
Mgll absorption is indeed the galaxy responsible for 
the absorption. 

We will demonstrate that galaxy disks with realis- 
tic aspect ratios contribute a geometric cross-section 
for Mgll absorption competitive with that of spheri- 
cal halos. We also examine the possibility that ob- 
served absorption properties may actually be gov- 
erned to a large extent by gas distributed in a disk 
geometry. In §3 we describe Monte-Carlo simula- 
tions of populations of galaxies with absorbing gas 
distributed in spherical and in disk geometries, and 
compare these models to the observed impact param- 
eters, luminosities, and equivalent widths in the SDP 
survey. A careful consideration of the selection pro- 
cedures, given in §4, is essential to the interpretation 
of our model results. In §5, we present simple pre- 
dictions, to help distinguish whether absorption often 
arises primarily in the disks themselves, based upon 
forth-coming high spatial resolution HST images of 
the absorbing galaxies and from high resolution spec- 
tra of the absorption lines. The basic results of our 
Monte-Carlo study of the cross-section of Mgll ab- 
sorption are summarized in the concluding §6. 

2. The Statistical Contribution of Disk Ge- 
ometries 

It is expected that some fraction of the Mgll 

absorbing galaxies likely result from lines of sight 
through the disks of Milky Way-like galaxies. High 
velocity clouds in the Halo have a covering factor 
of at most 38% [through half the halo (Savage, Lu, 
& Sembach 1995, Bowen, Blades, & Pettini 1995a)], 
and thus by themselves do not provide a near-unity 
covering factor. Further, Bowcn, Blades, & Pettini 
(1995a) find that a significant fraction of the high ve- 
locity clouds are kinematically consistent with disk 
co-rotation, so that they all may not really be a halo 
component per se, but instead be more closely asso- 



ciated with the disk. 

Since roughly 80% of field galaxies^ have disks, 
we consider the Mgll absorption cross-section con- 
tributed by a randomly oriented population of disks. 
In fact, a disk embedded in a sphere of the same ra- 
dius {Rd/Rfi = 1) presents a significant cross-section. 
An infinitely thin disk yields a lower limit of Y2 for 
the relative probability for intersecting a disk relative 
to a sphere. When the disk thickness is increased, 
or the condition Rd/Rt = 1 is relaxed, this proba- 
bility increases rapidly. We define the disk thickness 
T = h/Rd, where h is the half-thickness of the disk, 
and derive the relative probability for intersection of 
a population of randomly oriented disks with respect 
to spheres of radius R^, 

P{T,Rd/Rh) = {T + 0.5)(^^^ . (2) 

As illustrated in Fig. 1, to achieve equal cross-section 
contributions from a population of disks and spheres, 
the absorbing gas in a disk need only extend ^ 20% 
beyond that in a spherical structure for reasonable 
thickness. Values ofT = 0.2 — 0.3 are plausible for 
the outer, warped extensions of galaxy disks (Diplas 
& Savage 1991), since a disk with a scale-height of 
5 kpc would likely have a significant amount of ab- 
sorbing gas at twice that height. Furthermore, it is 
not unreasonable to conjecture that Mgll absorbing 
gas in these regions could be as extended or even more 
extended than a distribution of Mgll absorbing gas 
in a halo, given that the HI necessary for photoioniza- 
tion shielding is likely more concentrated in the plane 
defining the disk. 

It is also important to consider what the effective 
covering factor, down to some equivalent width limit, 
would be for a population of disks. We suggest that 
the equivalent width is likely to increase with disk in- 
clination because the increased pathlength through 
an inclined disk yields increased column densities 
and a larger velocity dispersion of absorbing clouds, 
since the number intercepted is proportional to path- 
length. A simple increase in column density increases 
Wo(A2796) most effectively when absorption falls on 
the linear part of the curve of growth. A curve 
of growth analysis reveals that, near the equivalent 
width limit of 0.3 A, absorption is on the linear part 
of the curve of growth for 6 > 20 km s~^. Doppler 

^The presence of Mgll absorption does select the population of 
field galaxies (Bergeron & Boisse 1991, Steidel 1993a). 
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parameters on this order are quite cornnion in the 
intermediate resolution (~ 30 km s~^) spectra of Pe- 
titjean & Bergeron (1990), and are consistent with 
the idea that low resolution data sample the velocity 
dispersion of a few or more individual clouds along 
the line of sight. For Wo(A2796) > 0.6 A, it is even 
more likely that the number of clouds intercepted and 
their velocity dispersion dominates the value of the 
equivalent width. Petitjean & Bergeron (1990) find 
that VFo(A2796) still scales linearly with the number 
of resolved subcomponents. 

Ultimately, the extent of the disk contribution to 
the overall cross-section for Mgll absorption depends 
upon the detailed physics of sharp HI edges (gas 
dynamical processes, ionization field and strength, 
galaxy mass and shape, external pressure) (Maloney 
1993, CorbelU & Salpeter 1993, Dove & ShuU 1994). 
Locally, extended HI disks (beyond the optical ra- 
dius) are common [cf. the recent review by Irwin 
(1995) and references therein], but only in a few 
cases has the Af(HI) distribution been mapped down 
to a sensitivity below the expected sharp edge at 
which the gas converts from mostly neutral to mostly 
ionized (Corbelli, Schneider, & Salpeter 1989, van 
Gorkom et al. 1993). In NGC 3198, the iV(HI) dis- 
tribution is observed to drop sharply with radius for 
iV(HI) < 5 X 10^9 cm-2 (van Gorkom et al. 1993), 
but is still 10^^ cm-2 35.40 j^-p^. ^pj^jg j-g^pj^j 

decline is likely to slow at larger radii. Bowen, 
Blades, & Pettini (1995b) find that there exists a 
correlation between the B magnitude of spiral galax- 
ies and their HI radii measured from 21-cm emis- 
sion. Their application of the Holmberg-like relation 
i?(HI) = R*{Lb/L*jjY to local HI disks provides evi- 
dence that (3 becomes less steep as the iV(HI) thresh- 
old is decreased, and approaches the values required 
to match the intermediate redshift Mgll absorbers, 
namely (3 ^ 0.2 [for B luminosity (S95)] at the Ly- 
man limit, iV(HI) 10^'^-^ cm~^. Although it may 
be an extreme example, we also note that in the M81 
interacting group there is a large covering factor of 
material with A^(HI) > 10^" cm~^ beyond a radius of 
40 kpc around M81 itself, which is spread throughout 
the group in a flattened distribution (Yun, Ho, & Lo 
1994). 

In the Milky Way, N{YLl) > 10^^ cm^^ j^^- 
tocentric radius of 30 kpc in the direction I = 130° 
(Diplas & Savage 1991). By extrapolating this value 
using their observed radial scale length of 5.5 kpc in 
this direction, we obtain iV(HI) = 10^^ cm~^ at a 



radius of ~ 60 kpc. HST spectra of various QSOs, 
Active Galactic Nuclei, and Halo stars, looking in 
various directions through the Galaxy, show Mgll 
absorption arising in the disk (Savage 1993, Lu, Sav- 
age, & Sembach 1994, Cardelli, Sembach, & Savage 
1995, Sembach, Savage, & Lu 1995), though the ra- 
dial extent of the Mgll itself is not observed directly. 
In some cases, the absorbing gas could actually be as- 
sociated with an A'^(HI) even smaller than 10^'^ cm~^ 
[at high redshift Mgll is found to be associated only 
with systems that are optically thick at the Lyman 
limit, but at low redshift the ionizing background is 
reduced and some Mgll systems are observed with 
tll somewhat less than unity (Bergeron et al. 1994)]. 
Thus, in the Milky Way, it might even be expected 
that the Mgll disk extends well beyond 50 kpc. To 
make all this plausible, metals would have to exist 
at these large radii, either by the action of galactic 
fountains, by infall of debris from satellites, or be- 
cause they were present in the material from which 
the galaxy was formed. 

Based upon these arguments, we suggest that the 
cross-section contribution of thick and/or extended 
disks and their related material (eg. co-rotating clouds) 
may well be large enough to account for the major- 
ity of Mgll absorbers. Realistically, absorbers are 
probably quite diverse, including contributions from 
LMC-Iike satellites, streaming tidal debris formed in 
merging events, occasional tll low surface brightness 
galaxies, and even early-type galaxies of sufficient K 
luminosity [the latter two having been observed in the 
SDP survey (Steidel, private communication)]. 

3. Modeling Absorbing Galaxies 

We perform Monte-Carlo simulations that allow us 

to predict the numbers and properties of absorbing 
and non-absorbing galaxies in a sample of randomly 
generated fields centered on QSOs. We have designed 
the simulations so that our results can be compared 
to observations. Since the SDP survey represents the 
most comprehensive observational program to date, 
we briefly describe their work [but also see Drinkwa- 
ter, Webster, and Thomas (1995)]. 

Drawing primarily from the Sargent, Steidel, & 
Bokscnbcrg (1988) and Steidel & Sargent (1992) sur- 
veys, SDP obtained optical and IR images of 51 QSO 
fields in which 58 Mgll absorbers had been detected 
(0.3 < z < 0.9) and of 25 "control" fields in which the 
QSO spectra exhibit no absorption to a completeness 
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limit of VFo(A2796) = 0.3 A. The images were taken 
over a large range of band passes and the detection 
levels were targeted to ~ 2 magnitudes fainter than 
present day L* for the known absorption redshifts. 
A combination of filters were used in order that the 
rest-frame B magnitude of the galaxies could be mea- 
sured. Additionally, K band images were obtained 
for all fields. Confirmation of any absorbing galax- 
ies identified in the imaging portion of the program 
required follow up spectroscopy. Unless a field had 
several identified candidates, spectra were obtained 
in a pattern from small angular separation to larger 
until a match was found for the redshift seen in ab- 
sorption. This follow-up spectroscopy in most fields 
is nearly complete within a 8 — 10" radius of the QSO. 
For control fields, redshifts were obtained for nearly 
all galaxies in the images within 10" of the QSO. Pre- 
liminary results from this study are summarized in 
S95. 

For our model fields, the redshifts, luminosities, 
and impact parameters are known for all galaxies 
within 10" of the field center. The fraction of absorb- 
ing and non-absorbing model galaxies is thus repre- 
sentative of an unbiased sample of QSO fields. Our 
procedure is as follows: The galaxy redshift is se- 
lected at random in the range 0.3 < z < 0.9, roughly 
matching the observed distribution. The redshift de- 
termines the physical size of the region (corresponding 
to a radius of 10") within which this galaxy will be 
located^. The galaxy luminosity Lk is selected from 
a Schechter function, 

P{LK)dLK ~ {LK/L*j,r exp [-{Lk/L*^)] (ILk, 

(3) 

with a = —1.0. This is derived from the value, 
a = —0.7, observed for the luminosity function of 
the absorption selected sample (S95), which is then 
weighted inversely by the absorbing cross-section area 
of a galaxy to obtain the true luminosity function. 
The galaxy is located randomly within the 10" region 
surrounding the QSO, giving the impact parameter 
D. The "observed" equivalent width Wq along the 
line of sight is then determined based on the assumed 
geometry and model parameters. If the generated 
Wo exceeds the 0.3 A observation threshold, then the 
galaxy is designated as an "absorber" , and otherwise 
as a "non-absorber" . 



'The conversion assumes go = 0. All sizes quoted throughout 
this paper would be ~ 80% smaller for go = 0.5. 



3.1. Spherical Models 

Invoking a spherical geometry, S95 suggested that 
the relatively sharp observed cut off in the size of 
Mgll absorbing regions may be provided by the de- 
crease of the pressure of diffuse gas that confines the 
clouds responsible for Mgll absorption within the 
halo. When the pressure decreases below a thresh- 
old value, which scales with the mass of the galaxy, 
the column density of HI in an individual cloud is 
no longer sufficient to shield Mgll and keep it in its 
singly ionized state. This effect roughly sets a mini- 
mum column density for individual clouds, for as the 
pressure decreases outward from the galaxy center, 
N{m) of the clouds decreases (Mo 1995). This col- 
umn density limit apparently corresponds to an equiv- 
alent width of 0.15 — 0.30 A, since weaker absorbers 
are rarely observed (Steidel & Sargent 1992, Churchill 
1996). To investigate this supposition, we ran a series 
of photoionization models, CLOUDY (Ferland 1988), 
and computed the equivalent widths of model clouds. 
These models were plane parallel slabs illuminated 
from one side with continuum slope Uox = -1.5. For 
the range of ionization parameters —4.0 < logF < 
— 1.0, the average equivalent width was 0.13 A for 
tll ^ 1 clouds with logA/'(HI) ~ 17.5 cm"^. 

Empirically, the distribution of equivalent widths 
is roughly consistent with a power law (Steidel & Sar- 
gent 1992), 

n{Wo)dWQ ~ Wo^dWo, (4) 

where S = 1.65 ± 0.90, and thus many absorption 

systems have total Wo near this limiting equivalent 
width. The relationship between the sum of the col- 
umn densities and the observed equivalent width de- 
pends upon the distribution of velocities of clouds 
along the line of sight. Intermediate resolution spec- 
troscopy reveals that the equivalent width increases 
almost linearly with the number of components (Pe- 
titjean & Bergeron 1990). Thus, it is reasonable to 
assume that for the low resolution spectra, the in- 
dividual clouds are sufficiently dispersed in velocity 
space so that their equivalent widths add to a total 
observed equivalent width Wq for the Mgll absorber. 

The challenge for designing realistic spherical cloud 
models is to produce a near-unity covering factor 
when a typical line of sight at large impact is likely 
to pass through only one or two clouds. If this is the 
case, then some lines of sight will certainly fail to pass 
through any clouds observable in Mgll, resulting in 
non-absorbers. 
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3.1.1. Design and Parameters 

Guided by this general picture, we design models of 
a distribution of clouds in a spherical halo, governed 
by the following parameters: (1) a mean number of 
clouds iVc(O) along a line of sight through the halo 
with impact parameter D = 0, (2) a power law depen- 
dence of the equivalent width of an individual cloud 
on the radius from sphere center, W{R) ^ i?^", that 
would result from a decrease in the individual cloud 
column densities due to decreasing pressure with in- 
creasing R, (3) a variable cloud size factor, given as 
the variance cr of a Gaussian distribution centered on 
W{R), and (4) a mean minimum limiting equivalent 
width VFmin for an individual cloud. 

For a line of sight at impact parameter D, the num- 
ber of clouds Nc{D) is determined by scahng A^c(O) 
by the pathlength through a sphere of radius R{Lk) 
[given by Eq. (1)]. In each model we choose either 
Wmin = 0.15 A or 0.30 A. For the former, the mean 
number of clouds at R{Lk) is set equal to two, and 
for the latter to one, so that the mean Wq at R{Lk) 
is 0.30 A in either case. The value of R{Lk) is deter- 
mined using /? = 0.15 to match observations (S95), 
however, we adjust the normalization R* to mini- 
mize the sum of the number of non-absorbers below 
and the number of absorbers above the R{Lk) line 
given by Eq. (1). This iterative process results in 
self-consistent models. The actual number of clouds 
along the particular line of sight is selected from a 
Poisson distribution with mean Nc{D). The loca- 
tions of the clouds are chosen at random along the 
line of sight (no dependence of cloud density on R), 
and the power law index a is used to compute W{R). 
The actual individual cloud equivalent width Wc is 
selected from a Gaussian distribution with variance a 
centered on W{R). The "observed" equivalent width 
Wo is determined by summing the Wc, an approxima- 
tion justified by the nearly one to one correlation seen 
between equivalent width and the number of subcom- 
ponents in an absorption line (Petit jean & Bergeron 
1990) . A successful model should recover the observed 
equivalent width distribution n{Wo)dWo. In order to 
roughly constrain our models, we assume a power law 
distribution and determine the least squares fit slope 
5 to the distribution of model Wq. 

3.1.2. Absorption Properties 

Table 1 lists the various series of model parameters 
and the fractions of absorbers above /above and below 



/below the R{Lk) boundary [shown in Fig. 2 and given 
by Eq. (1)] for an unbiased sample of fields. In fact, 
the fraction /below gives the covering factor within 
radius R{Lk)- From models SI S4, we see that 5 is 
strongly dependent on the value of a, but that a has 
little effect on the covering factor. There is a strong 
dependence of 5 on the mean number of clouds along a 
line of sight (models S5-S7). The dependence of these 
gross properties on a is minimal (compare models S3 
and S6). By increasing the number of clouds A"c(0), 
the covering factor /below can be increased such that 
it is quite close to unity, and a can be chosen so that 
the observed value oi 5 ^ 1.65 is obtained. In these 
spherical cloud models, large Wq result from multiple 
clouds along a line or sight and/or from a large value 
of a. Table 1 illustrates the balance between the pa- 
rameters iVc(O) and a. For a = 2, we can obtain 
(5 = 1.6 (enough large Wq) with a mean of only two 
clouds at D = (model Sll). For a ~ 0.5, a mean of 
three clouds is necessary (model S2). 

The covering factor could be further increased by 
using an even larger A'c(O), but then there would be 
too many large Wq values {5 would be too small for 
any a > 0). This problem is less severe if Wmin is 
reduced from 0.30 A to 0.15 A (compare models S7 
and SIO), but the lack of observed small Wmin values, 
discussed above, indicates that 0.15 A is really a limit 
on how far down this parameter can be pushed. For 
models S2, Sll, and S12, we set iVc(O) = 3, 2, and 
6, and tuned Wmin, oi and a to match the observed 
n(Wo)(iWo distribution. These three models will be 
used in the following sections to illustrate variations 
between different realizations of spherical cloud mod- 
els. 

The predicted distribution of impact parameter D 
versus K luminosity is illustrated in Figs. 2(a-c) for 
three different randomly selected realizations of model 
S2. To facilitate some level of comparison to Fig. 2 of 
S95, we plot 58 random absorbing galaxies. The num- 
ber of non-absorbers (63) presented on this figure is 
the mean number that would arise in the same num- 
ber of fields that yield a mean of 58 absorbers. At first 
glance, the lack of a clean demarcation of absorbers 
and non absorbers aromid the R{Lk) line appears 
inconsistent with the SDP data. As we will address 
in §4, the appearance and subsequent interpretation 
of a D versus Lk diagram may be quite sensitive to 
observational selection procedures. In fact, models 
S2 and S12 succeed fairly well in that they yield a 
large value of /below and a small value of /above- How- 
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ever, wc note that this is not the case for all spherical 
cloud models that we can design. For example, Ta- 
ble 1 shows that in model Sll a much larger fraction 
(0.34) of the galaxies below the R{Lk) line will be 
non-absorbers and a much larger fraction (0.22) of 
the galaxies above the R{Lk) line will be absorbers. 
In these models, absorbing galaxies can exist above 
the R{Lk) line because of the Poisson fluctuations 
in the number of clouds along the line of sight and 
the Gaussian spread around the equivalent width of 
a single cloud. From Table 1, we infer that the pa- 
rameter Nc{0) is the most critical in determining the 
value of /above- For iVc(O) — 2 the value of Nc{D) will 
decrease slowly with D and the probability of detect- 
ing one cloud can be significant for an impact greater 
than D = R{Lk). 

In Figs. 3(a-c), we also illustrate the distributions 
of Wo versus D for the three realizations presented 
in Figs. 2 (a c). These diagrams illustrate that our 
approach to modeling "dumpiness" within the ab- 
sorbers manifests in a Wq — D distribution qualita- 
tively consistent with the data presented in Fig. 3 of 
S95. The scatter in the diagram at a given D value 
is produced by the variation of the number of clouds 
along a line of sight, and the various locations (in 
galactocentric distance R) at which clouds are inter- 
cepted. 

3.2. Disk Models 

Expected properties of absorbing gas in galaxy 
disks should be guided by impressions of the thick, 
warped outer disk regions of nearby spirals, as we 
have discussed in §2. There are likely to be clumps 
and irregularities ("clouds") in the disk (Irwin 1995, 
and references therein). The number of clouds and 
their overall velocity dispersion is likely to increase 
with increased pathlength through an inclined disk. 
Based upon the conclusions of Bowen, Blades, & Pet- 
tini (1995b) that the extent of Mgll absorbing gas in 
disks is likely follow the Holmberg-like relation for B 
luminosity with f3 ^ 0.2, we assume that absorbing 
gas in disks roughly follows Eq. (1) for K luminosity 
with (3 = 0.15. 

3.2.1. Design and Parameters 

The parameters for our disk models are: (1) a 
power law dependence of W{R) on the radius in the 
disk, given by 



where i?o.3 is the radius at which a line of sight per- 
pendicular to the disk would yield Wq = 0.3 A on av- 
erage, (2) a "dumpiness" factor, given as the variance 
(T of a Gaussian distribution centered on W{R), and 
(3) a "cut-off radius" i?cut beyond which W{R) effec- 
tively vanishes, which could be somewhat larger than 
i?o.3 for the less severe ionization conditions likely at 
low redshifts (Bergeron et al. 1994). 

In order to determine Wq for a given line of sight, 
the disk is randomly oriented with inclination i [with 
probability P{i)di ~ sm.{i)di\. The line of sight is 
located at a random angle 6 from the major axis of the 
ellipse representing the inclined disk projected onto 
the plane of the sky. The radial position where the 
line of sight pierces the plane of the disk is then 



R = D 



cos^ e + 



smf 
cosi 



1/2 



(6) 



W(ii) = 0.3 (i?/i?o.3)"" A, 



(5) 



The equivalent width W{R) is given by Eq. (5) scaled 
by the increased pathlength factor sec(i) to roughly 
account for the combined effects of higher line of 
sight column density and velocity dispersion of the 
absorbing clouds. The "observed" Wq is then selected 
from a Gaussian distribution with variance a cen- 
tered on W{R). From physical principles, we expect 
a sharp cut-off at some radius in the Mgll absorb- 
ing gas, so we assign a cut-off radius in the range 
1.5 < Rcut/Ro.3 < 2 (except for one model with 
■Rcut = C!0, chosen for purposes of illustration). As 
with the spherical models, we normalize Eq. (1) with 
(3 = 0.15 by minimizing the sum of the number of 
non-absorbers below and the number of absorbers 
above the R{Lk) line. 

3.2.2. Absorption Properties 

In Table 2, we illustrate the effect of varying the 
disk model parameters. Note that for a popula- 
tion of infinite disks (model Dl), for which we fix 
R* = 38/i~^ kpc, there are no non-absorbing galaxies 
expected below the R{Lk) line, but c;learly disks must 
have some cut-off. As i?cut is reduced (models Dl- 
D3), the fraction of non-absorbing galaxies below the 
R{Lk) line increases, but the sensitivity in the regime 
-Rcut/-Ro.3 ~ 2 is small. From models D4-D6, we note 
that the fraction of non-absorbers below the R{Lk) 
line is not very sensitive to the choice of a. An inho- 
mogeneous population with 1.5 < Rcut/Ro.3 < 2 and 
1.0 < a < 2.0 would yield a similar fraction of non- 
absorbers /below ~ 0.70-0.75, provided that all galax- 
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ies can be characterized by R{Lk) = R*{Lk / L*j^ f'-^^^ . 
Again, a reahstic single population model must be 
consistent with the observed distribution n(Wo)rfWo 
with power 5 ~ 1.65. Models D7-D9 illustrate the 
level at which a, a, and i?cut can be adjusted while 
still successfully recovering the n(Wo)dWo distribu- 
tion and showing how the absorption properties vary 
between these models. We select these models for 
further discussion in the following sections. 

Using the same criterion as for spheres, we present 
the distributions of D versus Lk for three realiza- 
tions selected from model D8 (58 absorbers and 58 
non-absorbers) in Figs. 2(d-f). The fraction of galax- 
ies above the R{Lk) line that produce absorption 
is larger than /above = 0.15 for all cases in Table 
2. If i?cut is reduced further (below 1.5). /above can 
be reduced, but then the expected fraction of non- 
absorbing galaxies below the R{Lk) becomes much 
larger. It is difficult to design a realistic disk model 
that will not produce a non-negligible fraction of ab- 
sorbing galaxies above the R{Lk) line. Physically, 
this is because of lines of sight that pass through the 
disk at large impact parameters, but also at large in- 
clinations. 

We also present the Wq versus D distribution for 
58 absorbers from model D8 in Figs. 3(d-f). As with 
the spherical models, the Wq versus D distribution 
is qualitatively consistent with the S95 data, except 
that our disk models yield a deficiency of small Wo 
at small D. This is an artifact of our model design, 
which does account for increased pathlength through 
highly inclined disks using a geometric scaling, but 
does not account for material physically distributed 
above or below the disk plane (eg. thick disk or warp 
material). In a more realistic model, lines of sight at 
small impacts that graze highly inclined galaxies with 
some absorbing gas distributed at 5 — 10 kpc about 
the disk plane [as is likely for the Galaxy (Diplas & 
Savage 1991)] would also exhibit absorption, filling 
in the lower left corner of the — D distribution. 
It is interesting to note that in S95 most (if not all) 
of the absorbers with D < 15/i^^ kpc are Damped 
Lya (DLA) absorbers. This may imply that most 
DLA absorbers are highly inclined disks, and that, in 
view of our model interpretations, extended material 
above and below the disk is further implied by their 
detection. 



3.3. Ruminations on an Unbiased Sample 

The bottom line is that, regardless of the overall 
geometric cross-section, it is difficult to design mod- 
els of absorbing gas which (1) consist of physically 
distinct clouds as opposed to a smooth distribution 
of gas that would exhibit a tighter Wq — D relation- 
ship than is observed, (2) recover the distribution of 
equivalent widths n{Wo)dWQ, and (3) yield covering 
factors near unity down to the required equivalent 
width limit within some well defined luminosity de- 
pendent galactocentric radius. 

We note that none of our spherical cloud mod- 
els can yield few non-absorbing galaxies at impact 
parameters D < R{Lk), qualitatively reproduce the 
observed Wq — D distribution, and be made consis- 
tent with the observed equivalent width distribution. 
The basic problem is that gaps in covering result from 
instances of zero clouds along the line of sight. This 
could be adjusted by increasing the number of clouds, 
but then the equivalent width distribution is skewed 
toward large values. Decreasing the minimum equiva- 
lent width of an individual cloud could help with this 
discrepancy, but this would not be consistent with the 
small observed number of absorbers with Wo < 0.3 A. 
For disks, we note that taking into account the in- 
creased pathlength for inclined disks does in fact re- 
duce the number of non-absorbers at small impact 
parameters. Indeed, our disk absorber models yield 
similar fractions of non-absorbers below the R{Lk) 
line as do our spherical absorber models. In the S95 
sample, only three absorbing galaxies are observed 
above the R{Lk) line. Although some spherical cloud 
models appear to be relatively consistent with this re- 
sult, all of our disk models predict larger numbers of 
absorbers at large impact parameters. 

We must address the fact that our models yield 
no definitive boundary between absorbing and non- 
absorbing galaxies in the R{Lk) — Lk plane. It is 
important to take into account that the results we 
have discussed so far apply to an "unbiased sample" 
of QSO fields. For our unbiased sample, there does 
appear to be a difference between the predictions for 
spherical and disk models. It is in principle possi- 
ble to design a spherical cloud model that predicts 
very few absorbing galaxies above the R{Lk) line. 
Whether such a model is practically realistic, effec- 
tively whether there is a sharp cut-off in radius for the 
existence of absorbing clouds, depends upon the de- 
tails of ionization physics and on the variations in the 
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confining pressure and cloud propcrti(^s with radius. 
However, none of our disk models produce a negli- 
gible fraction of absorbing galaxies above the R{Lk) 
line, and we argue that all disk models would produce 
some absorbing galaxies at large impact parameters. 

In a unbiased sample of QSO fields (chosen with- 
out regard to the presence or non presence of absorp- 
tion) in which all galaxy redshifts and luminosities 
are measured to some limiting rest Lk, our "clumpy" 
models predict that there is no definitive boundary 
between absorbing and non-absorbing galaxies. We 
will show in §4, that a careful "application" of the 
SDP selection procedures to our model fields brings 
our results, for both spherical and disk models, into 
plausible agreement with observations. 

4. Confronting Observations 

The SDP sample is not an unbiased sample of QSO 
fields, but includes 51 "absorber fields" that exhibit 
Mgll in the redshift range 0.3 < z < 0.9 and 25 "con- 
trol fields", in which absorption is not detected. In 
order to better understand the effects of biasing in 
the sample, we studied Monte-Carlo fields of spheri- 
cal models S2, Sll, and S12 and of disk models D7, 
D8, and D9. An important (and not all that cer- 
tain) parameter is the mean number of galax;ies per 
field (within the redshift range). A'fioid can be es- 
timated from the number of absorbers observed per 
unit redshift, roughly dN/dz = 0.85 ± 0.2 (Steidel 
& Sargent 1992) in the range 0.3 < z < 0.9, which 
translates to iVios = 0.51 ±0.12 absorbing galaxies 
per line of sight. Prom Tables 1 and 2, we find that 
our simulations predict that the fraction of absorb- 
ing galaxies per fine of sight is 0.46 < /abs < 0.54. 
Thus, the predicted average number of galaxies per 
field is iVficid = 0.51//abs = 0.95 - 1.10. We adopt 
-Afield = 0.95, but note that the observational con- 
straint on A^ios allows some flexibility in this quantity. 

We generated 10000 simulated SDP surveys, each 
having 51 absorber fields and 25 control fields. The 
galaxies themselves were randomly selected from the 
models discussed in §3 (including redshifts). Based 
upon the presence or non-presence of absorption 
above the 0.3 A survey limit, we designated a field as 
an "absorber" or "control" field. Predictions based 
upon both the spherical and disk models are quite 
similar, and we thus discuss the simulation results 
more generally. For A^fieid = 0.95 we find that 38- 
41% of the fields (Fabs = 0.36 - 0.40) contain at 



least one absorbing galaxy. This is a much smaller 
fraction of absorbing fields than selected in the SDP 
sample (51/76 = 0.67), which may be telling us some- 
thing about our models or about a selection bias in 
the SDP survey. Simulated absorber fields contain 
on average 1.24-1.28 absorbers and 0.44-0.51 non- 
absorbers [c;ontrol fields contain the same number of 
non-absorbers (0.44-0.51), as they should since the 
presence of an absorber does not aff'ect the proba- 
bility of finding other galaxies in the field at differ- 
ent redshifts]. On average, the six geometric models 
yield 63-65 absorbers and 33-39 non-absorbers. The 
observed number of non absorbers (14) is not consis- 
tent (smaller) with the predictions of our models, and 
the number of absorbers observed (58) is somewhat 
improbable, on the low end of the tail of the distribu- 
tion. In Fig. 4(a), we illustrate the simulated distribu- 
tions of the number of absorbers and non-absorbers 
for model D8 in an "SDP survey" , in which all galax- 
ies within a 10" field have been included. In this 
figure, we also present the number of non-absorbers 
with impact parameters less than 8" which lie below 
the R{Lk) line. For this sub-sample, the predicted 
number non-absorbers approaches the observed num- 
ber, providing a powerful illustration of how sensitive 
observational results are to observational complete- 
ness and selection effects. 

4.1. Effects of Covering Factor and Selection 
Bias 

It turns out that the numbers of absorbers, non- 
absorbers, and the fraction of absorbing fields Fabs 
are sensitive to the parameter iVfieid- In Table 3, we 
illustrate the effect of decreasing A'fieid for models D8 
and S12 (other models produce similar numbers). As 
illustrated in Fig. 4(b) for model D8 and iVfieid = 0.65, 
the 51 model fields yield a mean of A^abs = 59 — 61 
absorbers, consistent with the 58 observed by SDP. 

A value of iVfieid = 0.65 is not implausible. If 
the observed dN/dz is about l-tr smaller, we recover 
A^iioid = 0.65. Alternatively, the estimate of TVfioid 
scales inversely with the fraction of galaxies in the 
fields that produce absorption, /abs- If we had de- 
signed models that yield nearly unit covering factor, 
we would have increased /abs and provided a means to 
reduce A/geid- However, this would yield an increased 
number of absorbing galaxies observed in 51 absorb- 
ing fields. Galaxy models that yield higher covering 
factors cannot be made consistent with the number 
of absorbing galaxies observed by SDP. Provided our 
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models statistically represent the absorption proper- 
ties of real galaxies, i^abs of an unbiased sample of 
fields may be as small as 0.28, and the SDP sample 
may have selected heavily in favor of absorbing fields. 

However, even with A^seid ~ 0.65, our models pre- 
dict that the observed number of 14 non-absorbing 
galaxies is somewhat improbable. Since the observa- 
tional procedures of SDP would often miss galaxies at 
large angular separations, it is more instructive to fo- 
cus on the fact that SDP find only two non-absorbing 
galaxies below the R{Lk) line separating absorbers 
and non-absorbers. We will show that this observa- 
tional result can be statistically reconciled with our 
simulations for both spherical and disk models once 
we account for (1) biasing in the choice of fields and 
in the practical procedure employed by SDP to ob- 
tain follow-up redshifts of galaxies, and (2) the un- 
certainty that even a small fraction of the absorbing 
galaxies may have been mis-identified (see §1). 

4.2. Reconciling Selection Effects 

We turn our attention to the fact that a fair num- 
ber of the non absorbing galaxies in the SDP sample 
do not have measured redshifts, and are therefore not 
represented in Fig. 2 of S95. In practice, the follow- 
up redshift measurements of galaxies identified in 
the SDP fields were implemented using the following 
practical considerations (Steidel, private communica- 
tion). In control fields, redshifts were obtained for all 
galaxies within the 10" fields. However, in absorber 
fields, non-absorbers with large impacts are prefer- 
entially missing from the sample, since the search 
was sometimes stopped at ~ 8" once a galaxy with 
the redshift seen in absorption had been found. In 
nearly all cases, redshifts were obtained for all galax- 
ies within 8". 

In the context of these procedures, we ask the rel- 
evant question - "for our models, what is the proba- 
bility of detecting only a few non-absorbing galaxies 
below a statistically determined R{Lk) line in a D 
verses Lk diagram?" . In order to address the; par- 
tial coverage in the 8" < £> < 10" zone, we bracket 
the SDP procedure with two "observational scenar- 
ios": (1) all non absorbers with 10" of the QSO are 
included, and (2) only those non-absorbers within 8" 
of the QSO are included. In Table 4, we present the 
predicted mean number of non-absorbers below the 
R{Lk) line based upon these two scenarios for the 
various models and for the A^fieu values 0.95 and 0.65. 
Prom the N ir and A^ // entries, note that there is 



not much sensitivity to the chosen observational sce- 
nario. Nonetheless, in order to best emulate the ob- 
servational procedures, we applied scenario (1) to the 
25 control fields and scenario (2) to the 51 absorber 
fields. We then computed P(< 2), the probabilities 
of observing two or fewer absorbers below the R{Lk) 
line. We find that both sphere and disk models give 
non-negligible probabilities for observing two or fewer 
non-absorbers in the SDP sample. There is a 31% 
chance (for spherical model SI 2) that SDP would de- 
tect two or fewer non-absorbers below the R{Lk) line 
and a 5% chance for disk model D8. 

Since the SDP observational procedure was to 
search outwards in the absorber field until a galaxy 
was located at the redshift seen in absorption, it is 
possible that some "absorbing" galaxies have actu- 
ally been mis-identified [though, for several reasons, 
it is believed that few mis-identifications have been 
made (S95)]. It is likely that for each mis identified 
galaxy, a non-absorbing galaxy below the R{Lk) line 
and an absorbing galaxy above the R{Lk) line would 
not be represented in the survey. Thus, in addition to 
P{< 2)), we have tabulated P{< 4) to examine the ef- 
fects of possible mis-identifications. The probability 
of observing four or fewer non-absorbers is not negli- 
gible for models S12 (71%), S2 (40%), and D8 (26%), 
for Afield = 0.65 (corresponding to dN/dz = 0.54). 
This is further illustrated by the distribution function 
for model D8 presented in Fig. 4. The probability of 
detecting four or fewer non-absorbers is not negligible 
in all three disk models, ranging from 16-26%. 

We briefiy consider the number of absorbing galax- 
ies above the R{Lk) line. In Table 5, we present the 
predicted numbers of absorbing galaxies above the 
R{Lk) line and the probabilities of observing fewer 
than four such galaxies in an SDP sample. Both ob- 
servational scenarios, eg. including all galaxies out to 
the field limit of 10" or including only those within 
8", are presented. In particular, model S12, which 
is characterized by a large Ac(0) and a small Wmin, 
results in a large P{< 4). As such, it yields the prop- 
erties of a "sharp boimdary". For disk models, there 
are an abundance of absorbers above the R{Lk) line 
in a sample of 10" fields. For Afidd — 0.65, model D8 
predicts a mean of 11 absorbing galaxies above the 
line, but only 5 of these are within 8" of the center of 
the field and thus certain to be included in the obser- 
vational sample. If the observational procedures have 
lead to a bias against detecting absorbing galaxies at 
8-10" then all models have a non-negligible P(< 4). 



Charlton & Churchill 



Mgll Absorbing Galaxies 12 



4.3. Discussion 

Based upon our models and a study of the cov- 
ering factor issues implied by the Holmberg-like re- 
lationship for absorbing gas cross-section, we cannot 
distinguish whether Mgll absorption arises primarily 
within clouds distributed in a spherical halo or dis- 
tributed in a flattened thick disk. Both of these geo- 
metric models appear plausible with the existing data 
once we include uncertainties due to possible mis- 
identifications and a possible bias in field selection 
of too few non absorber fields. The number of non- 
absorbing galaxies observed at small impact param- 
eters depends strongly on how far a sample of QSO 
fields is from an unbiased sample. We conclude that 
the small observed number of non-absorbers in the 
SDP sample does not necessarily imply unity cover- 
ing factor. Once observational biases and uncertain- 
ties are taken into account, the covering factor can be 
made consistent with 70-80%. Indeed, if the covering 
factor does approach unity, it is difficult to reconcile 
the small number of absorbing galaxies SDP observed 
in their 51 fields with models that match all other ob- 
servational constraints. 

It is clear to us that our models are idealizations 
of the universe of real galaxies. Our disk models 
do not take into account the fact that some num- 
ber of halo clouds are surely present around a re- 
alistic disk galaxy. As a result, disk models yield 
some non-absorbers for cases where the line of sight 
passes too far out in an inclined disk to produce ab- 
sorption (beyond i?cut), but some of these lines of 
sight might yield absorption in a galaxy with clouds 
in an extended halo. Likewise, our spherical "halo" 
models do not contribute a disk component, which 
would increase the number of absorbers. The impli- 
cation, if we are to obtain 58 absorbers in 51 fields, 
is that the covering factor of the halo proper in such 
model galaxies would be reduced below ^ 70%. Ad- 
ditionally, our models neglect the presence of satel- 
lite galaxies. Satellite galaxies, merging or otherwise, 
may contribute a significant cross-section for absorp- 
tion (cf. York et al. 1986, Wang 1993). From dy- 
namical considerations, it may be expected that ab- 
sorption properties may also correlate with the pres- 
ence of satellites, their distance from the primary, and 
the direction of their rotation with respect to the pri- 
mary (prograde/retrograde). Local spiral disk galax- 
ies are known to have an average of one to two satel- 
lites (Zaritsky et al. 1993). Interestingly, for their 
sample, the disk luminosity does not correlate with 



the number of satellites, so that absorption due to 
satellites would not be expected to correlate with 
the luminosity of the primary. However, they do see 
the "Holmberg-effect" , in which the projected satel- 
lites show some evidence for having an excess near 
the minor axis of the primary. Roughly, in the case 
of the Milky Way, the LMC has Mb = -18 (and 
Nh = lO^^cm"^ at i? = 7 kpc) so it could have a ra- 
dius for Mgll absorption of ~25 kpc (Bowen, Blades, 
& Pettini 1995b), thus covering as much as ^ 25% 
of the halo (not including the stream!). The effect of 
accounting for satellite absorbers in our models would 
yield an increase in predicted numbers of absorbing 
galaxies at all impact parameters. Again, the impli- 
cation, if we are to obtain 58 absorbers in 51 fields, 
is that the covering factor of the halo proper around 
galaxies would be reduced well below ~ 70%. This 
would represent a significant shift in the design of our 
models. 

It is highly unlikely that a single pre-dominant ge- 
ometric cross-section can be invoked to predict the 
presence or non-presence of absorption. For one, it 
is not clear to what degree merging history plays a 
role in geometrically distributing absorbing gas. Sim- 
ilar merging scenarios at different redshifts might 
be expected to produce different absorption prop- 
erties, given redshift evolution of the meta-galactic 
UV background and the influence of ionization on 
the strength of absorption and the dissipative prop- 
erties of the gas. Also, it is not yet clear at what 
level disk/halo interactions, such as galactic fountains 
(Bregman 1980) or chimneys (MacLow & McCray 
1988), contribute to the gas content of extended halos. 
Since the gas dissipation rate is density dependent, it 
is possible that the statistically observed cross-section 
of absorbers is a function of Mgll column density. 
Gas that dissipates on a few Gyr timescale would 
likely infall (arising from merging, extra-galactic ac- 
cretion, or disk/halo processes) onto the disk, con- 
tributing to its thickness and affecting its chemical 
evolution. Gas that fails to dissipate, yet maintains 
tll ~ 1 at a minimum, would likely remain dis- 
tributed in the halo, its absorption strength being 
governed by the galactocentric pressure. This com- 
ponent of the material may not always remain bound 
to the galaxy. This picture has the prediction that, 
as seen in high resolution spectra, low column clouds 
should exhibit the full range of velocities across the 
absorption profile, often detectable in the; wings of 
stronger lines, or suppressing the continuum between 
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strong lines, whereas high c;olunni c;lou(ls should have 
an internal velocity dispersion consistent with a pop- 
ulation of randomly oriented disks. The consequences 
for inferring the presence or non-presence of absorp- 
tion from a known galaxy with a given impact pa- 
rameter are not clear. It may be that the statistical 
geometric cross-section is sensitive to the equivalent 
width threshold and the epoch from which the sam- 
ple of galaxies is drawn. An Wo threshold of 0.3 A 
at 0.3 < z < 1.0 may sample a relative contribution 
from both halos and disks. 

5. Testing Absorber Cross— Section and Kine- 
matics 

There are simple predictable kinematic and ab- 
sorption properties that may be useful for discern- 
ing between these two geometries. Both tests require 
high spatial resolution imaging of the galaxies iden- 
tified as absorbers. The kinematic tests also require 
high resolution spectra of the absorption. Both sets 
of data will soon be available [(HST images) Steidel 
& Dickinson, private communication; (HIRES/Keck 
spectra) Churchill (1996)]. 

5.1. Disk Geometry and Associated Kine- 
matics 

Absorption associated with galaxy disks is likely to 
exhibit the characteristic signature of rotating galax- 
ies (cf. Lanzetta & Bowen 1992). As such, correlations 
between the strength and variation (spread) in the ob- 
served absorption with increased galaxy inclination 
provides a statistical test to discern absorber geome- 
try. In Fig. 5(a), we show a scatter plot of the dis- 
tribution of absorber equivalent widths Wq as a func- 
tion of galaxy inclination i. Note that nearly edge-on 
galaxies arc predicted to have a median Wq ~ 2-4 A. 
If absorption arises primarily in the disk, the absorp- 
tion exhibits a larger mean Wq and larger spread as 
the galaxy becomes more inclined. If absorption is 
primarily generated in the spherical halos of galax- 
ies, there is no such correlation predicted from our 
models. 

From HST images, it will also be possible to mea- 
sure the angle of the line of sight relative to the po- 
sition angle of the projected galactic major axis (the 
equivalent of our 9). In Fig. 5(b) we illustrate that, 
statistically, if absorption arises in a disk geometry 
there should be twice as many absorbers along the 
major axis as along the minor axis. Non-absorbers 



exhibit the opposite behavior. In a simplistic spher- 
ical case, there is no relationship between the clouds 
and the disk, and thus the relative number of ab- 
sorbers to non-absorbers exhibits no dependence on 
the position angle, 6. 

Establishing statistical kinematic trends will rely 
on accurate measurements of the absorption line sub- 
components and of the absorbing galaxy's systemic 
redshifts. For our disk models, we have computed 
the difference between the line of sight velocity seen 
in absorption and the galaxy systemic velocity. For 
complex absorption profiles with several subcompo- 
nents, one can think of the velocity seen in absorp- 
tion as an optical depth weighted mean (Hobbs 1973). 
We assume this weighted mean of the absorption ve- 
locity can, on average, be represented by the line of 
sight component of the TuUy Fisher velocity given 
by Lk/L*j, = [ViLK)/V{L*j,)]^-^ (Pierce & Tully 
1988, Forbes et al. 1995). The diff'erence between the 
absorption velocity and the systemic velocity of an 
Lk galaxy is then 

AV/V{Lk) ^ iD/R)cosesini. (7) 

We illustrate in Fig. 6(a) that large AV/V{Lk) can 
result when the galaxy has a relatively large inclina- 
tion. The distribution spread results from the various 
of the line of sight (eg. the small values at large i 
result when the line of sight passes close to the mi- 
nor axis). Indeed, as we show in Fig. 6(b), there is 
a strong statistical relationship between AV/V{Lk) 
and 9. In this case, the spread is a result of various 
galaxy inclinations. For establishing these trends ob- 
servationally, it is important to keep in mind that the 
accuracy of galaxy systemic redshifts can usually be 
measured with accuracy no greater than ^ 30 km s"-'^, 
but even with this uncertainty, the general trends are 
easily discernible for the expected ~ 200 km s~^ ro- 
tation speeds of galaxies. 

5.2. Spherical Geometry and Infall Kinemat- 
ics 

A very attractive picture, in which absorbing gas 
is distributed in more or less a spherical geometry, is 
radial infall of intergalactic material. Steidel & Sar- 
gent (1992) observed that redshift evolution of the 
co-moving cross-section of Mgll absorbers depends 
strongly upon the equivalent width limit of the sam- 
ple. The observed distribution of equivalent widths 
of Mgll absorbers changes such that the absorption 
strength becomes weaker with time. The argument 
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for infall (Stcidol 1993b) is based upon the conclu- 
sion that this evolution is not provided by changes in 
ionization conditions or metallicities, since Mgll ab- 
sorption almost always exhibits saturation. Since the 
sizes of the overall absorbing regions apparently do 
not evolve [for Wo(A2796) > 0.3 A, the cross-section 
is consistent with no evolution], the inference is that 
either the numbers or the velocities of the clouds are 
decreasing with time. In a spherical halo model, these 
ideas are further supported by the expected short life- 
times of T ~ 10^ K clouds moving in a virialized halo 
(Mo 1995), so that a reservoir for replenishment of 
gaseous material is required to account for the non- 
evolving halo size over a Hubble time. 

If absorption does arise from clouds distributed in 
a spherical halo with infall kinematics, the velocity 
width of the overall absorption profile should exhibit 
a trend with impact parameter. For smooth mass 
distribution the trend is quite simple, as we illustrate 
with the solid line in Fig. 7 for an L"^ galaxy with 
Rh = R*- The total velocity width is produced by 
the absorbing gas moving with the extreme line of 
sight velocities toward and away from the observer. 
This width is a maximum at £> = and decreases to 
no absorption for D > Rh- In the case of individ- 
ual infalling clouds, the profile velocity width is due 
to the maximum velocity difference between the two 
most kinematically extreme clouds. For the case of 
a constant infall velocity (no dependence upon R), 
we compute the profile velocity width [normalized by 
V{Lk)] for model S2, based upon our Monte-Carlo 
simulations. We illustrate the results as a scatter di- 
agram in Fig. 7. In cases of lines of sight that pass 
through only one absorbing cloud no point is plotted. 
Note that the general trend given by the smooth dis- 
tribution is present, though the scatter is quite severe. 

The features that are most relevant to future obser- 
vations are the deficit of small profile widths at small 
impact parameter and of large profile widths at large 
impact parameter. For the case of spherical infall, we 
re-cmphasize that siich observed trends would be cor- 
roborated if there was no evidence for a dependence 
on the galaxy orientation. However, it is not unlikely 
that infall would be influenced by the non-symmetric 
potential of the galaxy close in (small D), even if it 
were dominated by dark halos at large D. If so, then 
the spherical case presented in Fig. 7 would still ap- 
ply for large impact parameters, but the small impact 
parameter trends may exhibit additional scatter and 
show some connection to the galaxy orientation. 



6. Conclusions 

In this paper, we have been motivated by the cen- 
tral question: what is the predominant geometrical 
cross-section of galactic gas which is known to give 
rise to Mgll absorption in the spectra of QSOs? It 
has been our aim to ascertain if a single geometric 
cross section can be inferred as the general shape 
that governs the presence or non-presence of absorp- 
tion in a given galaxy of known properties interven- 
ing to a QSO. To address this question, we have per- 
formed Monte-Carlo simulations of Mgll absorption 
produced in either a spherical or in a disk geome- 
try. The models were designed to recover the observed 
statistical properties of Mgll absorbers, namely the 
distribution of Wo(A2796), given by n{Wo)dWo = 
W^'^-^^dWo (Steidel & Sargent 1992), the scatter in 
the Wo — D plane, and a near unity covering factor. 
We assume throughout that each model represents a 
single population of absorbers, and used the D — Lk 
plane and the Holmberg— like relation (Eq. 1) to com- 
pare our models to currently available observations. 
We find that: 

1. A population of randomly oriented disks can 
produce a geometric cross-section comparable 
to that of spherical halos. The relative con- 
tribution is a function of both the disk thick- 
ness and its extent beyond the halo proper (see 
Fig. 1). Additionally, lines of sight that impact 
highly inclined disks at large radii likely yield 
absorption with Wq > 0.3 A, due to the in- 
creased column density along the line of sight, 
the increased velocity dispersion of the clouds 
in the disk, and material that extends even as 
little as 5 kpc above and below the disk (such 
as with warps). 

2. Both spherical and disk models for which the 
parameters are tuned to recover n(Wo)(iIFo, 
yield similar fractions of non-absorbers at im- 
pact parameters smaller than the "boundary" 
given by R{Lk) = R*{Lk / L^^f-^"^ ■ Spherical 
models can be designed to yield small fractions 
of absorbers above the R{Lk) boundary, but 
disk models predict larger numbers due to lines 
of sight at large impact parameters through 
highly inclined disks. As a caveat, we note that 
the largest observed Wq may be sampling galaxy 
pairs (Churchill 1996), which implies that the 
distribution derived from single galaxies would 
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have a steeper slope. The observed trend for 
Mgll absorbers to become less strong with time 
(Steidel & Sargent 1992) may have some inter- 
esting implications for the questions we have at- 
tempted to address in this paper. One possi- 
ble view is that a "pre-dominant geometry" in 
which Mgll absorption arises could be chang- 
ing with time (evolving population), based upon 
the notion that numerous and kinematically di- 
verse optically thick clouds could be associated 
with halos or merging events at redshifts of 
z > 1 and could then be evolving into thick 
and warped disks at redshift z < 1 with rem- 
nants in the form of high velocity clouds still 
accreting at the present epoch. Such accretion 
is known to be occurring in the Galaxy [cf. the 
recently discovered dwarf galaxy in Sagittarius, 
Ibata, Irwin, & Gilmore (1994)], whose halo ap- 
pears to not be dynamically mixed (Majewski 
1996). Indeed, such mergers are predicted to 
leave detectable "streaming tails" or "moving 
groups" in the halo for timcscalcs of a Gyr or 
more (Johnston, Spcrgcl, & Hernquist 1995). 

3. Our models qualitatively reproduce the observed 
Wq — D distribution (Fig. 3 in S95), which sug- 
gests that our scheme to model the "clumpy 
structure" of Mgll absorbing gas in galaxies 
is appropriate, and provides some support that 
the covering factors predicted by our models are 
likely to be fairly accurate. We find that the 
small number of observed non-absorbers in the 
SDP sample does not necessarily imply near- 
unity covering factor to the equivalent width 
limit Wo(A2796) = 0.3 A, but can be consis- 
tent with covering factors as small as 70-80%. 
In the case of disks we can think of an effec- 
tive covering factor determined by considering a 
population of randomly oriented disks. If Mgll 
absorbing galaxies do indeed have covering fac- 
tors approaching those predicted by our mod- 
els, then the 58 absorbers detected by SDP in 
51 QSO fields imply that the number of galax- 
ies per 10" QSO field may be as small as ~ 0.65 
in the redshift regime 0.3 < z < 0.9 in which 
case dN/dz is ~ 1-a smaller than the observed 
value. Spherical cloud models with larger cover- 
ing factors cannot be made consistent with the 
observed n{WQ)dWo and Wq — D distributions 
and also reproduce as few as 58 absorbers in 51 
fields in our "SDP survey" simulations. 



4. Since our models arc cither entirely "halo" or 
"disk" , they both likely under predict the num- 
ber of observed absorbers. The implication is 
that the halo "component" of a model incorpo- 
rating both disk and halo would necessarily have 
a covering factor < 70%, if we are to obtain 58 
absorbers in our 51 simulated absorber fields. 
In our simulations, we find a gap in the Wq — D 
plane at small Wq for highly inclined disks at 
small impact parameters [D < 15/i~^ kpc). We 
attribute this to the fact that our model disks 
are not physically extended above or below their 
plane. Since most or all observed absorbers in 
this impact regime are DLAs (S95), and since 
we infer that these absorbing galaxies indeed 
have "thick" disks based upon our models, we 
suggest that many of the small impact DLAs 
arise from highly inclined thick disks. 

5. Our models, which are an unbiased sample of 
QSO fields in which all galaxy redshifts are 
"measured" to a limiting Lk, generally yield 
no clear boundary between absorbing and non- 
absorbing galaxies. This is consistently true 
for disk models, but for spherical models the 
"distinctness" of the boundary depends strongly 
on parameter choices. A paucity of absorbers 
above the R{Lk) line could be used to argue 
against a disk geometry. But we note that a 
relatively larger number of absorbers above the 
line could be consistent with either geometry. 

6. Using the the observed number of non-absorbing 
galaxies that fall below the R{Lk) line in S95 as 
a discriminant between a population of spheri- 
cal (halo) absorbers and disk absorbers, we find 
that, once we account for (1) possible selection 
biases toward absorbing fields, and, (2) possible 
mis identifications of a few galaxies, we cannot 
definitively distinguish which geometry governs 
the presence or non-presence of absorption. We 
find that the probability for detecting only a few 
non-absorbers below the R{Lk) line is not neg- 
ligible for cither our sphere or disk models. The 
final results of the SDP survey, of which S95 is 
only a preliminary report, are likely to provide 
a demanding test for the validity our models. 
In principle, the number of absorbing galaxies 
that fall above the R{Lk) line may be used as a 
discriminant as well, with small numbers incon- 
sistent with a disk cross-section but consistent 
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with certain spherical cloud models. A com- 
plete unbiased sample would allow a more ro- 
bust measure of the fraction of absorbing galax- 
ies both above and below the R{Lk) line. 

7. Upcoming high spatial resolution images (HST) 
of the absorbing galaxies and high resolution 
spectra (few km s~^) promise to provide the 
necessary data to settle the issue of how Mgll 
absorbing gas is distributed in and around galax- 
ies. In particular, if absorption does arise in the 
disks of galaxies, both the disk inclination and 
the angle subtended between the location of the 
line of sight on the plane of the sky and the 
galaxy major axis are expected to be strongly 
correlated with presence of absorption and with 
the velocity difference between the absorption 
centroid and the galaxy itself (see Fig. 5). If 
absorption arises purely within halo clouds (in- 
falling or otherwise kinematically distinct from 
the disk), no such trends with disk orientation 
are expected. No doubt, there will be examples 
of each, and of both together in a single sys- 
tem. The real issue we have attempted to ad- 
dress is which one statistically governs the pre- 
diction of absorption or non-absorption, since 
this should yield important clues about the ori- 
gin and fate of gas for theories of galaxy evo- 
lution. For clouds infalling into spherical halos, 
the width of the overall absorption profile is ex- 
pected to decrease with increasing impact pa- 
rameter (see Fig. 7), though the scatter can be 
quite large in the regime 20 < f < 40 Un- 
til these observations and tests are investigated, 
we submit that the standard picture of a spheri- 
cal distribution of halo clouds around galaxies is 
not necessarily required by existing data. Since, 
at intermediate redshifts, Mgll absorption may 
be present in the disks of galaxies out to radii 
at which HI drops below 10^^ cm~^, we suggest 
that Mgll disks are quite likely to make a sub- 
stantial contribution to the population of Mgll 
absorbers. 

In the real universe it is likely that both spherical 
and disk geometries contribute to the Mgll cross- 
section, and that no single geometric configuration 
of discrete absorbing clouds can identified as being 
solely responsible for Mgll absorbers. If this is true, 
then the apparent boundary in the D — Lk plane, 
which provides compelling evidence suggestive of a 



mechanism for a K luminosity (mass) dependent cut- 
off in absorbing galaxies, may indeed result from small 
number statistics. 

As an example of how varied absorption conditions 
may be, consider three possible lines of sight through 
a disk galaxy: 1) A line of sight close (< 10-20 kpc) 
to the center of a modestly inclined spiral galaxy must 
pass through the disk and may pass through halo 
clouds or satellite galaxies. If the inclination is large, 
then for absorption to be strong, the disk would likely 
have extended material above or below the plane. 2) 
If the line of sight passes at larger impact parame- 
ter through a modestly inclined disk galaxy, then the 
question of the origin of Mgll absorption is more con- 
troversial. It becomes a question of whether the HI 
disk extends far enough beyond the optical radius of a 
galaxy and/or whether the halo has a significant cov- 
ering factor of clouds, high velocity or otherwise. The 
relative importance of these two contributions may 
change with redshift. 3) Very large impact lines of 
sight will pass too far out in an inclined disk galaxy for 
absorption from the disk itself to dominate. In these 
cases, an absorption line may still result from halo 
clouds or from satellites or satellite debris and would 
likely be optically thin. Though there are contribu- 
tions to absorption profiles from multiple components 
in galaxies and their surroundings, a general geomet- 
ric cross-section for absorption may still be definable. 
Chances are that this will be a function of the Mgll 
equivalent width threshold, or actually, of the column 
density "contour". Progress on defining the contri- 
butions of halos and disks to these cross-sections is 
forthcoming through high-resolution imaging, while 
the specific details of what galaxy components com- 
bine to produce absorption will rely on their signa- 
tures in high-resolution spectra. 
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Table 1 

Effect of the Variation of Parameters for Sphere Models 



Model Parameters Absorption Properties 



Model 


a 


a 




[A] 


R* [kpc] 


/abs 


/below 


/above 


<5 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


SI 


0.0 


0.2 


3 


0.30 


34 


0.46 


0.76 


0.08 


2.3 


S2 


0.5 


0.2 


3 


0.30 


35 


0.48 


0.78 


0.09 


1.8 


S3 


1.0 


0.2 


3 


0.30 


35 


0.47 


0.79 


0.08 


1.4 


S4 


2.0 


0.2 


3 


0.30 


36 


0.48 


0.80 


0.10 


1.2 


S5 


1.0 


0.4 


2 


0.30 


37 


0.46 


0.66 


0.22 


2.1 


S6 


1.0 


0.4 


3 


0.30 


35 


0.47 


0.78 


0.09 


1.4 


S7 


1.0 


0.4 


6 


0.30 


36 


0.53 


0.93 


0.03 


0.4 


S8 


0.0 


0.4 


6 


0.15 


37 


0.48 


0.83 


0.04 


3.5 


S9 


0.5 


0.4 


6 


0.15 


37 


0.48 


0.84 


0.04 


2.1 


SIO 


1.0 


0.4 


6 


0.15 


37 


0.49 


0.85 


0.04 


1.4 


SU 


2.0 


0.4 


2 


0.30 


39 


0.47 


0.66 


0.22 


1.7 


S12 


0.8 


0.4 


6 


0.15 


37 


0.49 


0.85 


0.04 


1.6 



Note. — The model designation is given in column 1. Column 2 lists the index for W{R) ~ for individual 

clouds. The variance a around this mean equivalent width W{R) is listed in column 3. Column 4 gives the mean 
number of clouds along & D = Q line of sight. The minimum equivalent width Wmin of an individual cloud is listed 
in column 5. Column 6 lists the value of R{Wo = 0.3 A) for galaxy. Columns 7, 8, and 9 give the fraction of 
all galaxies in unbiased fields that produce absorption, the fraction of galaxies below the boundary that produce 
absorption, and the fraction of galaxies above the boundary that produce absorption, respectively. In column 10 is 
the best fit slope of an assumed power law to the Wq distribution. 
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Table 2 

Effect of the Variation of Parameters for Disk Models 



Model Parameters Absorption Properties 



Model 


a 


a 




R* [kpc] 


/abs 


/below 


/above 


(5 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Dl 


1.0 


0.0 


oo 


38 


0.91 


1.00 


0.79 


1.5 


D2 


1.0 


0.0 


2.0 


34 


0.52 


0.75 


0.23 


1.9 


D3 


1.0 


0.0 


1.5 


37 


0.47 


0.70 


0.18 


1.8 


D4 


0.0 


0.0 


2.0 


27 


0.50 


0.69 


0.26 


3.0 


D5 


1.0 


0.0 


2.0 


34 


0.52 


0.75 


0.23 


1.9 


D6 


2.0 


0.0 


2.0 


35 


0.49 


0.75 


0.16 


1.5 


D7 


1.0 


0.4 


2.0 


38 


0.54 


0.71 


0.33 


1.7 


D8 


1.3 


0.0 


2.0 


34 


0.50 


0.74 


0.20 


1.8 


D9 


1.1 


0.2 


1.5 


40 


0.50 


0.71 


0.24 


1.8 



Note. — The model designation is given in column 1. Column 2 lists the index for 
W(K) ~ for individual clouds. The variance a around this mean equivalent width 

W{K) is listed in column 3. Column 4 gives the cut-off radius in units of i?(Wo = 0.3 A). 
The value of R{W^ = 0.3 A) for an L*j^ galaxy is listed in column 5. Columns 6, 7, and 8 
give the fraction of all galaxies in unbiased fields that produce absorption, the fraction of 
galaxies below the boundary that produce absorption, and the fraction of galaxies above 
the boundary that produce absorption, respectively. In column 9 is the best fit slope of an 
assumed power law to the Wq distribution. 



Charlton & Churchill 



Mgll Absorbing Galaxies 



Table 3 

Effect of the Variation of Density of Absorbing Galaxies 



Full Sample of Fields'' Absorber Fields'^ Control Fields'^ 

-Afield -^abs -^abs ^non ^ahs ^non -^abs ^non 



Model D8 



0.95 


0.38 


64 


36 


1.26 


0.48 





0.48 


0.80 


0.33 


62 


30 


1.21 


0.40 





0.40 


0.65 


0.28 


60 


25 


1.17 


0.33 





0.33 



Model S12 



0.95 


0.37 


64 


37 


1.25 


0.49 





0.49 


0.80 


0.32 


62 


31 


1.21 


0.41 





0.41 


0.65 


0.27 


59 


25 


1.17 


0.33 





0.33 



Note. — A^fieid is the input mean number of of galaxies per field. The fraction of fields 
that have at least one absorber is given by Fabs- For the full sample, the absorber fields, 
and the control fields, the A^abs and A^non denote the mean number of absorbers and non- 
absorbers, respectively. 

^Total mean numbers in a sample of 51 absorber fields and 25 control fields 
''Mean numbers in individual absorber fields 
■^Mean numbers in individual non-absorber fields 
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Table 4 

Expected Numbers of Non-absorbers "Below the R{Lk) Line" 



TVfieid = 0.95 ATfieid = 0.65 



Model 


^10" 


^8" 


P(< 2) 


P(< 4) 


^10" 


^8" 


P(< 2) 


P(< 4) 


S2 


9 


7 


0.018 


0.123 


6 


5 


0.109 


0.402 


Sll 


13 


12 


0.000 


0.007 


9 


8 


0.009 


0.080 


S12 


6 


5 


0.110 


0.412 


4 


3 


0.314 


0.712 


D7 


12 


10 


0.002 


0.022 


8 


7 


0.030 


0.170 


D8 


10 


9 


0.006 


0.052 


7 


6 


0.051 


0.260 


D9 


12 


10 


0.002 


0.023 


8 


7 


0.028 


0.163 



Note. — ^iq" is the mean number of non-absorbing galaxies below the line given by 
Eq. (1) in the text. N^n is the mean number of non-absorbing galaxies below the line that 
are within 8" of the QSO line of sight. P{< 2) is the probability of observing two or fewer 
non-absorbing galaxies below the line, and P{< 4) is the probability of observing four or 
fewer non-absorbing galaxies below the R{Lk) line. To mimic the selection procedures of 
the SDP sample, the probabilities are computed for the number of galaxies within 8" of the 
QSO line of sight in the 51 absorber fields plus the number within 10" for the 25 control 
fields. 



Table 5 

Expected Numbers of Absorbers "Above the R{Lk) Line" 



Model 




-Afield 


= 0.95 






^fiold 


= 0.65 






P(<4) 


^8" 


P(<4) 




P(<4) 




P(<4) 


S2 


6 


0.322 


3 


0.778 


5 


0.375 


3 


0.822 


Sll 


13 


0.001 


6 


0.233 


13 


0.003 


6 


0.279 


S12 


2 


0.934 


1 


0.992 


2 


0.953 


1 


0.994 


D7 


18 


0.000 


7 


0.157 


16 


0.000 


7 


0.199 


D8 


11 


0.008 


5 


0.410 


11 


0.016 


5 


0.464 


D9 


14 


0.001 


6 


0.238 


13 


0.002 


6 


0.289 



Note. — ^^q" is the mean number of absorbing galaxies above the line given by Eq. (1) 
in the text. N^n is the mean number of absorbing galaxies above the line that are within 
8" of the QSO line of sight. P(< 4) is the probability of observing four or fewer absorbing 
galaxies above the R{Lk) line. 
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Fig. 1 . — The effect of disk tliickness and radius on tlie geometric cross-section of a population of randomly oriented 
disks. The vertical axis is the relative probability of intercepting a disk with thickness T = h/Rd and radius Rd 
relative to that for intercepting a sphere of radius Rh- Here, h is the half-height of the disk. The four curves are 
labeled by their corresponding T. 
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Fig. 2. — The distribution of impact parameter versus normalized K luminosity for three random realizations of 
the spherical model S2 (a-c) and of the disk model D8 (d-f). Open circles denote absorbing galaxies and filled 
triangles denote non-absorbing galaxies. These realizations sample the number of fields required to produce 58 
absorbing galaxies, as observed by SDP. The model then yields the number of non-absorbers for this same number 
of fields, unbiased with respect to whether absorption occurs in the field or not. The solid line is the S95 best-fit 
to the SDP data, D = ^Sh'^lK / L^^f '^^ kpc, which we call the "R(Lk) line". 
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Fig. 3. — The distribution of "observed" Mgll equivalent widths Wo versus impact D for the models presented in 
Figure 2. — (a^c) Three randomly selected realizations of the spherical cloud model S2. In our spheres models, 
the scatter in Wo arises from clouds intercepted at different galactocentric distances and due to variations in the 
numbers of these clouds. — (d-f) Three randomly selected realizations of the disk model D8. For disks models, 
scatter is dominated by their random inclinations. 
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Fig. 4. — The distribution of the numbers of absorbers and non-absorbers predicted from model D8 in a sample of 
51 absorber fields and 25 control fields (the simulated SDP survey). From right to left, each histogram represents 
(i) the number of absorbing galaxies, (ii) the number of non-absorbing galaxies, and (iii) the non-absorbing 
galaxies with impact parameters less than 8" which lie below the R{Lk) line. To match observations, the mean of 
these distributions should be consistent with the observed numbers: 58, 14, and 2, respectively. — (a) Predicted 
distributions for mean of 0.95 galaxies/field. — (b) Predicted distributions for mean of 0.65 galaxies/field. The 
smaller number of galaxies/field provides a closer match to the observed numbers, although the numbers of non- 
absorbers are still on the low-end tails of the distributions predicted by our model. 
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Fig. 5. — Statistical correlations between absorption strength and galaxy inclination and orientation with respect 
to the line of sight are predicted if the absorption is in fact associated with the galaxy disk. — (a) The trend of 
increasing absorption strength Wq with disk inclination for model D8. The spread at fixed inclination is a result of 
a range of impact parameters. — (b) The predicted relative numbers of absorbing and non-absorbing galaxies as 
a function of the angle 9 subtended between the sky location of the line of sight and the major axis of the oriented 
galaxy disk. The absorbers are denoted by solid circles and the non-absorbers by open circles. Disk models yield 
absorbers preferentially at angles close to the major axis, and non-absorbers to the minor axis. Spherical cloud 
models, such as the ones we have designed in which the velocities and numbers of clouds have no connection to the 
galaxy disk, produce no correlations with disk geometry. 



1 



Charlton & Churchill 



Mgll Absorbing Galaxies 





30 60 90 

Inclination [degree] 



90 180 270 

9 [degree] 



360 



Fig. 6. — Statistical correlations between absorption kinematics and galaxy systemic velocity are predicted if the 
absorption arise in the extended disk. To observe such correlations will require accurate redshifts in absorption and 
for the host galaxy. — (a) For model D8, the weighted mean line of sight velocity of the absorption line relative 
to galaxy center verses disk inclination. The radial velocity is expressed in units of the rotation velocity of an Lk 
galaxy. The spread in observed velocity differences increases with disk inclination. — (b) The radial velocity of 
the absorption line relative to the center of the galaxy verses the angle 9 of the LOS relative to the major axis for 
model D8. Scatter at a given 9 results from the distribution of disk inclinations. 
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Fig. 7. — For a picture in which clouds are infalling into spherical halos with constant infall velocity, the profile 
velocity width (high resolution spectra) should be the maximum velocity difli"erence between the two kinematically 
extreme clouds. The maximum cloud velocities intercepted along a line of sight will be correlated with impact 
parameter as shown here. The individual points are from Monte-Carlo simulations of model S2. At very small 
impact parameter the majority of the points are at nearly the maximum value of 2V(Lk), since the motion of 
the clouds is directed along the line of sight. However, occasionally due to small number statistics all clouds are 
moving in the same direction so that the maximum velocity difference is near zero. The solid line illustrates a 
smooth distribution of infalling absorbing gas for an galaxy with Rh = R* ■ 
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